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Abstract—Plasma confinement in the field of a levitated dipole offers many advantages for magnetic fusion. 
MHD stability is obtained from compressibility which utilizes the large flux tube expansion of a dipole field. 
Such a device could be high beta, steady state, and exhibit good confinement properties. The large flux expan- 
sion will ease the difficulty of the divertor design. The configuration is ideal for electron cyclotron heating and 
controlled convective flow patterns may provide a mechanism for fueling and ash removal. 


1. INTRODUCTION 


The dipole magnetic field is the simplest and most 
common magnetic field configuration in the universe. It 
is the magnetic far-field of a single, circular current 
loop, and it represents the dominate structure of the 
middle magnetospheres of magnetized planets and neu- 
tron stars. The use of a dipole magnetic field generated 
by a levitated ring to confine a hot plasma for fusion 
power generation was first considered by Akira Haseg- 
awa after participating in the Voyager 2 encounter with 
Uranus [1]. Hasegawa recognized that the inward diffu- 
sion and adiabatic heating that accompanied strong 
magnetic and electric fluctuations in planetary mag- 
netospheres represented a fundamental property of 
strongly magnetized plasmas not yet observed in labo- 
ratory fusion experiments. For example, it is well- 
known that global fluctuations excited in laboratory 
fusion plasmas result in rapid plasma and energy loss. 
In contrast, large-scale fluctuations induced by sudden 
compressions of the geomagnetic cavity (due to 
enhancements in solar wind pressure) or by unsteady 
convections occurring during magnetic substorms 
energize and populate the energetic electrons trapped in 
the Earth’s magnetosphere [2]. The fluctuations induce 
inward particle diffusion from the magnetospheric 
boundary even when the central plasma density greatly 
exceeds the density at the edge. Hasegawa postulated 
that if a hot plasma having pressure profiles similar to 
those observed in nature could be confined by a labora- 
tory dipole magnetic field, this plasma might also be 
immune to anomalous (outward) transport of plasma 
energy and particles. 


The dipole reactor concept is based on the idea of 
generating pressure profiles near marginal stability for 
low-frequency magnetic and electrostatic fluctuations. 
From ideal MHD, marginal stability results when the 
pressure profile, p satisfies the adiabaticity condition, 
(p V”) = 0, where V is the flux tube volume and y= 5/3. 
From gyro kinetics, marginal stability results when 
dF(u, J, y)/dyw = 0 where F(u, J, y) is the particle dis- 


tribution function, u is the adiabatic invariant, u = 
€; /2B, J is the parallel invariant, J = fvi dl and the flux, 


y is the third adiabatic invariant. The frequencies that 
correspond to these invariants are respectively Q, the 
cyclotron frequency, œ©,, the bounce frequency and , 
the curvature driven precessional drift frequency. Both 
of these conditions lead to dipole pressure profiles that 
scale with radius as 7°? while the adiabatic distribu- 
tion function, F(u, J) also implies a density dependence 
of n, « r^ and a temperature dependence T œ r®3, 


By levitating the dipole magnet in order to prevent 
end losses, conceptual reactor studies supported the 
possibility of a dipole fusion reactor [3, 4]. The dipole 
reactor concept is a radical departure from the better 
known toroidal-based magnetic fusion reactor con- 
cepts. For example, the most difficult problems for a 
tokamak reactor are the divertor heat dissipation, dis- 
ruptions, steady state operation, and an inherently low 
beta limit. Furthermore, the tokamak is subject to neo- 
classical effects and micro-turbulence driven transport. 
The dipole concept provides a approach to fusion 
which solves these problems. 


1. Divertor problem: The difficulty in spreading the 
heat load at the divertor plate is generic to concepts in 
which the magnetic flux is trapped within the coil sys- 
tem. By having the plasma outside of the confining coil 
the magnetic flux can be sufficiently expanded to sub- 
stantially reduce divertor heat loads. 


2. Major disruptions: A tokamak has a large amount 
of energy stored in the plasma current. The dipole 
plasma carries only diamagnetic current and is inher- 
ently free of disruptions. Experiments observe relax- 
ation oscillations, reminiscent of tokamak ELMs, 
which relax the pressure gradients when they become 
too large. 


3. Steady state: A tokamak is a pulsed device and 
current drive schemes that are required for steady oper- 
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ation appear to be costly. The dipole plasma is inher- 
ently steady state. 


4. Beta limit: Tokamak stability depends on the 
poloidal field which is less than the toroidal field by 
B,/Br ~ 1/qA with q the edge safety factor and A the 
aspect ratio. For a dipole there is a critical pressure gra- 
dient that can be supported (due to the compressibility 
term in the interchange stability criterion) and for a suf- 
ficiently gentle pressure gradient the dipole plasma 
resides in an absolute energy well and is stable up to 
local beta values in excess of unity. 


5. Transport and neoclassical effects: The trapping 
of particles in regions of bad curvature makes the toka- 
mak susceptible to drift frequency range trapped parti- 
cle driven turbulent transport. A dipole could, theoreti- 
cally, show classical transport. In addition a tokamak 
has a “neoclassical” degradation of transport that 
derives from the drifts of particles off of the flux sur- 
faces. In a dipole the drifts are toroidal and they define 
the flux surfaces. 


The chief drawback of the dipole approach is the 
need for a levitated superconducting ring internal to the 
plasma. Although this provides a challenge to the engi- 
neering of the device recent advances in high tempera- 
ture superconductors coupled with an innovative design 
concept of Dawson [5, 6] on the maintenance of an 
internal superconducting ring in the vicinity of a fusion 
plasma leads to the conclusion that this issue is techno- 
logically solvable. 


Utilization of electron cyclotron resonance heating 
(ECRH) and pellet injection in a levitated dipole per- 
mits a unique heating approach for the creation of 
fusion grade plasmas. The combination of closed field 
lines and high beta makes the dipole an ideal geometry 
for creating a high beta, hot electron plasma by the 
application of ECRH. Such plasmas have been demon- 
strated in open field line mirror machines as well as in 
levitrons and in the non-levitated dipole (or so-called 
“terrella’”) CTX experiment [7—10]. The density of the 
hot electron plasma is observed to be limited by micro- 
wave accessibility and the energy by relativistic detun- 
ing. Unlike ECRH mirror applications, the hot elec- 
trons created in a levitated dipole can only be lost by 
cross-field transport and if transport is close to classical 
the resulting loss rate will be small. The energy stored 
in the hot electron plasma can be transferred into a 
dense, thermal hydrogenic plasma by injecting pellets 
into the ECR heated plasma. 


We conclude that the continued advances in the 
understanding of magnetospheric plasmas coupled 
with the development of the technology of gyrotrons 
and of pellet injectors provides the fusion community 
with a new and unique approach to creating fusion 
grade laboratory plasmas. In addition the development 
of high temperature superconductors and of new con- 
ceptual approaches to the cooling of an internal levi- 
tated superconducting coil provide a potential path for 
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a substantially improved approach to magnetic fusion 
energy production. 


2. DIPOLE CONFIGURATION 


Figure 1 which shows the magnetic flux contours 
and the mod-B surfaces for a laboratory-size levitated 
dipole facility [11]. The levitated ring has a major 
radius R, and a minor radius a. The simplest vacuum 
chamber geometry would be a spherical tank and the 
last closed flux surface will be assumed to have a radius 
R„ Typically R,, 2 5R, to permit a substantial reduction 
of plasma pressure. The plasma pressure will attain a 
peak value at a midplane radius of Rọ. Between the ring 
surface (located at R = R, + a) and the peak pressure 
location there is good curvature and since the field is 
relatively high near the ring, cross-field transport is be 
expected to be low. Thus we expect high pressure gra- 
dients in this region. The schematic midplane pressure 
profile is shown in Fig. 2. 
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Fig. 1. Magnetic field and flux surfaces for a levitated dipole 
experiment. 
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Fig. 2. Midplane pressure profile shown schematically. 


Outside of the peak pressure region the pressure will 
decay in a region of “bad” curvature and MHD stability 
requires a relatively gentle pressure gradient. This 
results in a large volume of low field and of low density 
plasma. Thus a dipole confinement device would look 
like a small ring surrounded by a hot plasma that is cen- 
tered within a relatively large spherical vacuum cham- 
ber. Beta, the ratio of plasma to magnetic field pressure 
is a measure of the efficiency of the utilization of the 
magnetic field and since the field and pressure fall off 
together, B only decays slowly with radius. 


When a small uniform field is added in the direction 
of the dipole axis a separatrix will form which leads to 
two possible geometries: (1) When the uniform field 
adds to the outer dipole field the separatrix exhibits two 
symmetric nulls on the dipole axis as exhibited in 
Fig. 3. (2) When the uniform field subtracts from the 
outer dipole field the separatrix exhibits a field null on 
the outer midplane, as seen in Fig. 4. Within the last 
closed flux surface we expect to satisfy the MHD 
requirement (p V”) > 1. Outside the last closed flux sur- 
face the plasma will flow along field lines into the end 
wall. 


Beyond the separatrix heat will be conducted along 
the field lines and flow toward the end plates or limiter. 
A strong fanning of the diverted field lines is expected 
to greatly reduce the end wall power density. The 
plasma at the divertor plates is expected to be 1 to 
10 eV and at the symmetry point (opposite the end 
plate) 20 to 100 eV. As with a tokamak divertor the 
scrape-off layer plasma can be lead into a baffled cham- 
ber to permit pumping at high neutral pressure with a 
limited back streaming of the neutralized outflow. 


One straightforward way to heat a laboratory dipole 
experiment would be to utilize electron cyclotron reso- 
nance heating (ECRH). For an ECRH heated experi- 
ment the heating will take place at both the fundamen- 
tal and the second harmonic resonance field locations 
and the hot electrons will be localized between the flux 
surfaces that are tangent to the fundamental and the 
second harmonic mod-B surfaces. A typical laboratory 
embodiment of a levitated dipole could utilize 2.45 to 
28 GC ECR heating. The coil design used in Fig. [1] is 
constrained by the current density limit for a supercon- 
ductor, taken here to be 1 x 108 A/m?, and by the 
requirement that the resonant flux tube clears the coil. 


KESNER, MAUEL 


The vacuum field of a circular loop provides a useful 
analytic approximation. The vector potential for the 
field of a single loop of radius R, in cylindrical coordi- 
nates has the form 


_ wR 1/2 3 ` ; 
Ag = e(o) [d -0.5k )K(k)- E(k )], (1) 


with B = V x A and the magnetic flux surfaces are given 
by rAg. K and E are complete elliptic integrals of the 
first and second kind. The ring current is J and 
4R 
= — 
(R.+7r) +z 


For r > R, the field has a dipole dependence 
HoM 


3 
Tr 


B(r, 0) = “(1 + 3sin°0) (2) 
An additional vertical field term, Ag, = 0.5B,r term 


must be added to the vector potential (Eq. 1) to produce 
a diverted geometry. 


3. SCRAPE-OFF LAYER CONSTRAINT 


Since the pressure gradient must not exceed a criti- 
cal value the peak pressure, po, is determined by the 
scrape-off layer pressure: 


Po/ Psot = (Tsoi/To) ~ 10°10". 

Assuming that a fraction fp of power leaving the 
plasma is radiated we can balance the power loss from 
the plasma with the flow in the scrape-off layer: 


C1 - fr) PoVo 


= 2Pso1Âs01C5» 
VE 


with c, the scrape-off layer sound speed, Vo the effec- 
tive volume of hot plasma defined by pyVo = Í p dV, Aol 


the scrape-off layer cross-section area and Po (Pso1) the 
respective core (scrape-off layer) pressures. The zero 
subscript denotes values at the location of the pressure 
peak. Assuming that the scrape-off layer width is a 
finite number of ion gyro radii, we take A o1 = TYR soi) so1 
with P,o the scrape-off layer ion gyro radius and y the 
ratio of the scrape-off layer width to the ion gyro radius 


(y2 1). Taking Vo = TR, and imposing the critical tem- 
perature profile T œ R33 and y= 1 we obtain 


to" = 0.501 f ) ByRo (Rey (3) 
j l i NTioTeo Ro 


in MKS units with temperatures in eV. Consider two 
examples: 

1. A dipole with Rọ = 0.4 m, R,,, = 3 m, By = 0.4 T, 

fr = 9.5. For a hot electron plasma with T, = 300 KeV, 

PLASMA PHYSICS REPORTS Vol. 23 


No.9 1997 


PLASMA CONFINEMENT IN A LEVITATED MAGNETIC DIPOLE 


300 


400 
R,cm 


0 100 200 300 


Fig. 3. Configuration with separatrix null points in dipole 
axis. 


and Tọ = 200 eV, we obtain Tg ~ 0.1 sec. For a density 
of n, = 5 x 10! m” this corresponds to B, = 0.38 and the 
power requirement of 75 KW. This estimate of power 
ignores the stabilization of the hot electron interchange 
mode by the thermal electron component [12] which 
would permit a higher ratio of peak to scrape-off layer 
pressure and thereby reduce the power requirements. 
Therefore this heating power estimate may be consid- 
ered to be an upper bound. 


2. For Ry = 0.4 m, R.o = 3 m, To = Ta = 5 KeV, By = 
0.4 T, f= 0.5 we obtain Tg = 0.15 sec. For a density of 
n, = 2 x 10! this corresponds to B = 0.5 and the power 
required to sustain this plasma is 65 KW. 


For Tg < ge” from (3) the scrape-off layer can 
expand to accommodate a higher heat flow, i.e. y > 1. 


For Tg > we the heating power must be decreased, i.e. 
the temperature 7) must decrease for the heating power 
to balance the scrape-off layer power outflow with the 
divertor width fixed at one ion gyro radius. This behav- 
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Fig. 4. Configuration with separatrix null on outer mid- 
plane. 


ior results from the constraint imposed by the adiabatic 
profiles of pressure and temperature. 


We have assumed that T,,/T) = (Ro/Ryo))*®?. More 
generally the temperature profile is determined by the 
heating and thermal transport profiles. The relationship 
of the temperature to the density profile is also con- 
strained by stability considerations as discussed below. 
Furthermore the scrape-off layer temperature is deter- 
mined by the cooling processes that take place within 
the scrape-off layer including radiation and secondary 
electron emission. As the plasma flows along the 
scrape-off layer field lines toward the end plates it will 
cool and it can become collisional. The scaling changes 
when the mean free path within the scrape-off layer 
becomes less than the connection length and the heat 
flow along the field lines becomes predominantly con- 
ductive. For a collisional plasma the temperature will 
vary along the field line in the scrape-off layer (at con- 
stant pressure) due to parallel thermal conductivity. 
Additionally if the neutral pressure becomes suffi- 
ciently large (in the vicinity of the end walls) the pres- 
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sure can fall due to momentum exchange with neutral 
particles. 


For a collisional scrape-off layer the power balance 
becomes 


(1 — fr) Po Vo ~ PsotKil V sot 

Te joes 
and we note that for classical conductivity % ~ T”. 
Since collisionality retards the outflow of heat along 
open field lines the scrape-off layer temperature can be 
higher without degrading core confinement as com- 
pared with a collisional scrape-off layer. 


4. MHD AND ELECTROSTATIC MODE 
STABILITY 


The dipole reactor concept is based on the idea of 
generating pressure profiles near marginal stability for 
low-frequency MHD and electrostatic fluctuations. 
From ideal MHD, marginal stability results when the 
pressure profile, p(y) satisfies the adiabaticity condi- 
tion, 6(pV") = 0, where V is the specific flux tube vol- 
ume and y= 5/3 [13]. Hasegawa pointed out that when 
the invariants u and J are conserved, marginal stability 
of both MHD and drift modes results when F(u, J, 
y)/dy = 0, where F is the particle distribution function 
and y is the third adiabatic invariant. Both of these con- 
ditions lead to dipole pressure profiles that scale with 
radius as p « 729, 


At high beta MHD ballooning is expected to provide 
a stability limit. The stability of ballooning modes has 
been studied for both magnetospheric and fusion dipole 
applications [14-15] and these studies indicate a beta 
limit in excess of B > 1. For application to planetary 
systems or for supported dipole experiments one must 
take account of anisotropic pressure and also the 
boundary conditions where the field lines enter the 
poles. For fusion applications one is interested in a lev- 
itated dipole which has closed field lines and therefore 
we expect the pressure to be isotropic. For an isotropic 
pressure we can apply the standard MHD formalism. 


On the other hand, since confinement in a levitated 
dipole results from cross field transport, (and not pitch- 
angle scatter) it may be expected that, on a transport 
time scale, the distribution function will become isotro- 
pic. To lowest order the distribution function would be 
approximated by Fy = Fo(€, y) with € = (vi + vi )/2, 
the particle energy. It still remains true that low fre- 
quency unstable modes (œ < œ, Q.) would lead non- 
linearly to the adiabatic distribution function, dF(, J, 
w)/ow = 0. The subsequent collisional relaxation of the 
distribution function will preserve n, T, and n (N = 
dlnT/dlnn) and this has important consequences for 
micro stability. 
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Interchange Modes from Kinetic Theory 


When Fo = Fie, y) the kinetic analysis of electro- 
static drift modes proceeds from a traditional approach 
[16] beginning with the drift kinetic equation. A kinetic 
stability analysis of interchange modes produces simi- 
lar results to the MHD interchange criterion. In [16] it 
is shown that, interchange modes become stable when 
the following criterion is satisfied: 


(4) 
where @s = b x k, - V'Fo/mQ,Fo, is the diamagnetic 


drift frequency, @ =k, - b x (nv, b - Vb + LVB)/mQ, 
is the processional drift frequency, and b = B/|B|. In 
this analysis we have assumed œ < œ, with œ, the 
bounce frequency defined by œ, = v,d/ds and we have 
bounce averaged the drift kinetic equation. The bounce 
averaged quantities appear with an overline, & = 


[av FG; -B,0,) > 0, 


Ja dl/ vj. This yields the approximate stability criterion 


@, > (1 + N)@x which indicates that (at low p) stability 


requires that the local pressure scale length exceeds the 
radius of curvature. 


Trapped Particle Modes 


Trapped particle modes are low frequency electro- 
static modes (© < @, < Q,) that can localize in the 
outer part of the torus. For a tokamak @,/Q;. = 1/A < 1 
and as a result we can ignore resonances with the pre- 
cessional drift, @,, and obtain drift frequency fluctua- 
tions of order œ ~ @x. With a dipole we have the oppo- 


site ordering, i.e. ®,/@ > 1 and stability of interchange 
modes require @, > (1 + N)@x. The dipole ordering 


leads to the result that instability may be driven by a 
resonance of the wave with the processional drift 
motion, i.e. O ~ Wy. 

An analysis of the dissipative trapped ion mode that 
assumes collisional electrons and collisionless ions 
indicates a stability criterion 0, > N0/2 [16]. Like- 
wise an analysis of a collisionless trapped particle 
mode indicates a stability criterion @, > 7 [16]. 


Therefore we can conclude that, in a levitated dipole, a 
pressure profile that is stable to interchange modes may 
also be stable to trapped particle modes. 


N; Modes 


In a tokamak the so-called n; or “mixing” modes are 
electrostatic modes that couple ion acoustic and drift 
fluctuations. They are low frequency modes © < Q. 
but include the dynamics of the bounce frequency time 
scale. A dipole does not have magnetic shear and so we 
will refer to the early derivations of the n; mode which 
do not include magnetic shear. Antonsen et al. [17-18] 
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derives stability criteria from both a fluid and a colli- 
sionless kinetic theory. The simplest fluid theory indi- 
cates that instability occurs when n; > 2/3. The adia- 
batic distribution is characterized n; = VInT,/Vinn = 
2/3. The more detailed fluid derivation and the colli- 
sionless kinetic approach indicate that instability 
occurs when n; > 1. Thus adiabatic profiles which are 
characterized by N; = 2/3 are expected to be stable to 
these modes. 


Drift Cyclotron Modes 


Drift cyclotron modes are high frequency unstable 
modes (@ ~ Q) modes that are driven by temperature 
and density gradients. Pastukhov and Sokolov have 
evaluated the transport from these modes in the good 
curvature region near the surface of the levitated dipole 
[19, 20]. They show that the resulting transport would 
be severely limited by particle recycling at the surface 
of the internal coil. Because the surface of dipole is 
completely surrounded by a dense plasma, the net par- 
ticle flux to the ring must vanish. A cool, high-density 
sheath forms at the dipole surface which transforms the 
thermal flux into bremsstrahlung radiation. 


Stability of Convective Cells 


Experiments in multipoles have indicated that con- 
vective cells [21] can provide the dominant source of 
cross-field transport in shear-free systems. It is under- 
stood theoretically [22] that zero frequency convective 
cells are closely related to interchange modes and they 
will grow in regions of bad magnetic well, i.e. where 
(p V”) < 0. In addition it has been shown that convec- 
tive cells will exist in regions of good curvature when 
the heating is non-uniform [23]. In the Wisconsin octo- 
pole experiments convective cells were observed in 
regions of both good and bad curvature [22] regions 
and the convective cells were observed to decay. In 
these experiments the initial plasma was non-uniformly 
distributed within the chamber. These experiments also 
indicated that a small amount of shear caused a rapid 
decay of the convective cells. Additionally it was 
observed that small field errors can cause convective 
flow patterns in a shear-free configuration [24]. 


These results lead one to suspect that uniform heat- 
ing will be required in a levitated dipole in order to 
avoid the excitation of convective cells. Alternatively 
one can consider the addition of a small current on the 
dipole axis (or appropriate current drive) to create some 
shear and thereby eliminate these modes. 


On the other hand, convective cells can be usefill in 
a reactor to purge the ash. At marginal stability, pV’ = 
const, the exchange of flux tubes will transport particles 
but not energy. Additionally when n, « 1/V M; = 2/3) 
the exchange of flux tubes will not change the density 
profiles. Therefore for an adiabatic profile the con- 
trolled application of asymmetric heating can provide a 
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means for the adiabatic convection of the fusion prod- 
ucts. One might then consider the possibility that, in a 
reactor, a fraction of the synchrotron radiation that is 
emitted from the plasma can be reflected back so as to 
be non-uniformly reabsorbed and therefore drive con- 
vective cells that fuel and cleanse the plasma. 


Conclusions on Stability 


In general the radial profiles of heating and fueling 
determine the respective temperature and density pro- 
files. We can imagine heating and fueling profiles that 
generate marginally stable pressure profiles (p « V”) 
while maintaining n,; < 2/3. If we exceed the critical 
pressure gradient we would expect the plasma to 
expand unstably so as to broaden the pressure profile. If 
we exceed a critical value of n; we would expect an 
onset of micro-turbulence drive transport. This is in 
direct analogy to the operation of a tokamak. With pres- 
sure gradients and 1, appropriately bounded a levitated 
dipole may exhibit classical transport. 


5. ECR HEATING AND PELLET INJECTION 


Electron cyclotron resonance heating provides a 
straight forward method to build up a substantial stored 
energy in a hot electron plasma. This capability results 
from the predicted stability at high beta and the closed 
field line geometry of a levitated dipole. The energy 
stored in the hot electron plasma can be transferred into 
a high density, thermal hydrogenic plasma by injecting 
pellets into the ECR heated hot electron plasma. The 
energy transfer into the hydrogenic and the thermal 
electron species depends on a competition between re- 
thermalization and radial thermal transport. If transport 
is classical or near classical essentially all of the stored 
energy can be transferred into the hydrogenic and the 
thermal electron species. In this optimistic scenario 
D?He ignition could be obtained in an inexpensive 
manner in a relatively small reactor. 


If we consider, for example, a dipole plasma heated 
by a 28 GC gyrotron at the 1 T resonance surface (about 
100-200 KW of power) we may expect to produce a 
hot electron plasma with n, ~ 2-5 x 10'8 mœ at 200- 
500 KeV (0.2 < B < 1). The ability to transfer the energy 
depends on the competition of cross field transport and 
thermalization of the hot electron energy. For classical 
confinement, a large fraction of the energy could be 
transferred and a n, ~ 2 x 10!8 m, 200 KeV electron 
plasma would yield a 10 KeV hydrogenic plasma at 
ne ~ 2 x 10? m”. 


To study this process we utilize a zero dimensional 
model assuming classical re-thermalization between 
the suprathermal and the thermal species. Thermal 
plasma losses due to cross field transport characterized 
by Tg and Tp We solve the following simultaneous rate 
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Fig. 5. Solution to rate equations for electron and ion tem- 
peratures with Tg = a | for a= 0.1 m, and initial values of 
nep = 5X 108 m”, ng =5 x 109 m”. 


equations for hot electrons, T,,,, thermal electrons, T,, 
thermal ions, T;, and density, n, (assuming n, = n;): 


dT on = T= Toi + T; = Ten 


dt Toi Tin 
dT, =. Ten = Te + T; = T: T; 
dt 7 Tei Tei Tre 
dT; _ Ten = T; T; = Te T; 
dt Tn Tei Te; 
dn; n; 
dt T; 


with T, the classical rethermalization time, T;, the hot 
electron-ion rethermalization time, T,(T,;,.) the ion 
(electron) energy confinement and qT, the particle con- 
finement time. The energy confinement time is related 
to a mean transport coefficient by Tg = 67/4. 


For classical confinement essentially all of the 
energy can be transferred into the thermal plasma. One 
example would be to consider classical confinement 
with T; = a?/%_ for a = 0.1 m and initial values of n,, = 
5 x 10!8 nr, ng = 5 x 10° m>. As shown in Fig. 5 T, 
rises up to 12 KeV and 7; reaches 8 KeV. A second 
example, shown in Fig. 6 shows the solution to rate 
equations for with Tg = 20 ms T, = 200 ms and a = 
0.1 m. Even for substantially reduced ion confinement 
significant plasma parameters would be attained due to 
energy transfer from the relativistic electrons. Thus, by 
combining ECRH and pellet injection a modest dipole 
experiment might be capable of investigating the prop- 
erties of high temperature, high density, and high beta 
plasmas. 
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Fig. 6. Solution to rate equations for electron and ion tem- 
peratures with Tg = 20 ms, T, = 200 ms anda = 0.1 m. 


6. THERMAL EQUILIBRIA OF AN IGNITED 
D?HE PLASMA 


An estimate of the equilibrium profiles in an ignited 
D?He plasma requires a solution of a heat transport 
equation. This is a two-dimensional problem since the 
magnetic field strength varies along a field line (the 
density and temperature are flux functions) In a levi- 
tated dipole geometry a major portion of the volume of 
confined plasma is located near the outer midplane. We 
can therefore estimate confinement using a one dimen- 
sional model that uses fields evaluated at the outer mid- 
plane of the dipole. We assume D?He power production 
with a 20% D, 80% 3He mixture and furthermore 
assume that the dominant loss channels are 
Bremsstrahlung radiation and classical transport. The 
magnetic field also enters the Jacobian to describe the 
flux tube expansion and for classical transport kK, œ 
1/B?. We can write the thermal transport equation in 
flux coordinates using the magnetic scalar potential, x, 
as the angular coordinate, i.e. B = Vx. The flux tube 
average of the transport equation then gives: 


SI faxr’ 1 5 
(5) 
= [nony OV) ppe E fus(T) J Provan TI A 


On the right hand side (OV) ie is the DĉHe rate coef- 


ficient, Eps = 18.3 MeV and P,,em is the Bremsstrahlung 
radiation power given in [25]: 


Prrem = 5X 10” Tne [Zg +1.55 x 10°T, 


+7.15 x 10°T2) + 0.071(ZZ;n,/n,)T, 
+ 0.004147, ] 
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in W/m? with T in KeV and n, in m>. On the dipole 
midplane R = R, and (4) becomes: 


1 ð al 


A eee 
RoBoð Ro © Ro 


B |ƏR, = 


= [nony CO V) py Eus (T) gè Pron DISE. 


We have used k, œ 1/B? i.e. we have assumed classical 
transport. 


The flux average integrals, can be approximated by 
f (dl/B) R? œ RÉ and Í (dl/B) « R4. The scaling of 
(5) indicates that the temperature half-width, A; scales 


as Ar% ,/K,/H(T) with H(T) the sum of the heating 


terms from the right hand side of (5) Therefore a larger 
thermal transport leads to a broader temperature profile 
If we compare the total heat flux toward the ring vs 
toward the scrape-off-layer we can show that when 
Ar < Ror with Ro; the location of the temperature peak. 
The ratio of heat to the ring vs heat to the scrape-off- 
layer is Hp/H,,, = 1 + 6A7/Ror 


sol 

To solve Eq. (6) a density dependence must be 
assigned but since all terms go as n* this dependence 
only has a weak effect on the solution. (The pressure 
gradient, however, enters into the stability criteria and 
the pressure gradient broadens when the density profile 
broadens). The solution of (5) is a strong function the 
amplitude and radial dependence of x, which we have 
taken to be classical. Thus at a fixed beta, as the mag- 
netic field increases the peak fusion power density rises 
as B? but additionally the hot plasma region narrows 
due to reduced transport. 


7. DISCUSSION 


A dipole fusion reactor would consist of a single 
levitated circular magnet within a large vacuum cham- 
ber. The hot plasma core would encircle the levitated 
dipole coil forming a toroidal annulus. A large expan- 
sion region of cooler plasma extends outward from the 
dipole where the plasma pressure decreases with 
radius, R, approximately as R characteristic of the 
marginally stable profiles found in magnetospheres. 
Although the overall dimensions of the dipole fusion 
reactor may be large, the size of superconducting 
dipole magnet is small. Indeed, in the dipole reactor 
conceptual designs [3-5] the volume of the hot plasma 
core (40 m°) exceeds the volume of the levitated ring 
(20 m°). This feature of the dipole reactor (i.e. a larger 
plasma volume than the volume of the high-technology 
superconducting magnets, shield and structure) is in 
sharp contrast to the tokamak where the volume of the 
plasma is usually less than the volume of the surround- 
ing fusion island. (For example, in ITER, the plasma 
volume is 2500 m? and the volume of the magnets, 
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shield and structure exceeds 5000 m°). The dipole reac- 
tor concept also differs from the spherator [26] since 
the plasma profiles of the spherator are steep (i.e. they 
cannot be made stationary) and low-frequency fluctua- 
tions or convection cells may significantly degrade con- 
finement. 


Conceptual dipole reactor designs have been 
reported [3—4], and the use a dipole fusion reactor for 
space propulsion has been proposed [5]. In each of 
these designs, D—*He fuel was used instead of the more 
highly reactive D-T fuel in order to reduce the neutron 
flux to the levitated coil. Also reflectors were used to 
reduce synchrotron losses from the high-pressure and 
lower B plasma on the inside of the levitated dipole. 
The high B capability of the dipole reactor makes pos- 
sible the use of advanced and possibly aneutronic fuels, 
but the high temperatures required to burn these fuels 
necessitate steps to reduce synchrotron emission losses. 
The designs reported in [3] and [4] described compact 
and relatively low-power dipole reactors with large 
plasma expansion regions. A 20 MA dipole coil of 
radius of 1.8 m confined a plasma with peak B ~ 3 and 
generated 100 MW of fusion power A higher field, 
40 MA dipole with a denser plasma at the same B could 
generate 1000 MW. The plasma is heated to ignition 
with direct heating of the plasma core (using, for exam- 
ple, neutral beam injection). In [5], a much larger 
dipole reactor containing a 54 MA dipole having a 
radius of 6 m and producing 2000 MW of fusion power 
was considered for rocket propulsion The plasma 
expansion region was not as large and a relatively hot 
plasma was diverted to an annular gas-neutralizer to 
generate thrust In both [4] and [5], thermoelectric con- 
verters were located within the levitated dipole, and 
they provided the power to drive refrigerators for the 
superconducting magnets. Designs of the supercon- 
ducting magnets and shields in [3] and [5] illustrate the 
feasibility of reactor-sized dipole magnets using 
present-day multi-filamentary Nb3Sn conductors. 


Closely related studies of plasma confinement in 
multipole fields have been a subject of significant 
research efforts [27, 28]. Recent theoretical work 
includes the study of the equilibrium of two dimen- 
tional internal coil configurations [27-30]. 


Finally, a supported dipole laboratory experiment at 
Columbia University has studied the detailed phase- 
space evolution of dipole-trapped energetic plasma in 
the presence of intense drift-resonant fluctuations [9- 
10]. In this experiment, an “artificial radiation belt” (a 
population of 5-50 KeV energetic electrons) is pro- 
duced with microwave heating which excites hot-elec- 
tron interchange instabilities when the pressure gradi- 
ent sufficiently exceeds the marginally stable profiles 
envisioned for the dipole reactor. Fluctuations leading 
to global chaotic transport as well as thin, localized 
regions of stochastic drifts are observed, and these 
observations have been used to verify models of colli- 
sionless radial transport in dipole magnetic fields. 
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